Generating Gravitational Waves After Inflation 
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I review two mechanisms by which gravitational waves can be generated at the end of inflation: preheating, 
and gravitons Hawking radiated during the decay of very small primordial black holes. These mechanisms are 
contrasted with the gravitational waves during inflation, and may provide a window into the physical processes 
that govern the end of the inflationary phase. 



1. INTRODUCTION 

It has long been known that an inflationary 
phase in the very early universe sources a back- 
ground of gravitational waves. The amplitude of 
this background is tightly constrained by Cosmic 
Microwave Background [CMB] data [314] placing 
nontrivial constraints on the inflationary param- 
eter space [S]. The temperature anisotropy of 
the CMB is measured with much higher accuracy 
than the polarization anisotropy (although con- 
straints on the latter are rapidly improving f6J) 
but ultimately polarization measurements will 
provide the strongest constraints on the primor- 
dial tensor background [7]. The amplitude of 
this signal is strongly correlated with the physical 
scale of inflation, and there is no guarantee that 
inflation occurs at a high enough scale for this pri- 
mordial signal to be observable, even in an ideal 
experiment. In the longer term, space-based in- 
terferometers might make direct detection of this 
background at wavelengths of thousands or even 
millions of kilometers, but these comoving scales 
still left the horizon a number of e-folds before 
inflation ended. 

Recently, substantial attention has been paid to 
physical processes at the end of inflation which 
generate gravitational waves - either preheat- 
ing [1|8|9|10|11|12|13] , or Hawking radiation from 
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primordial black holes [2]. In contrast to the 
usual inflationary signal, these backgrounds are 
sharply peaked, and have frequencies correspond- 
ing to "laboratory" (or, at their lowest, solar 
system) scales, as opposed to the astrophysical 
scales associated with the CMB. Their amplitude 
can be substantial, but the technical challenges 
associated with the detection of high frequency 
gravitational waves suggest that it may be some 
time before these backgrounds are directly con- 
strained. Conversely, these backgrounds are sen- 
sitive to physical processes occurring at the end 
of inflation, and the coupling between the infla- 
ton field(s) and the rest of the matter sector, and 
could provide a unique window into the primor- 
dial universe. 



2. GRAVITATIONAL WAVES FROM 
PREHEATING 

Parametric resonance and preheating follow- 
ing inflation has been studied in great detail: 
see, for example |14ll5ll6ll7ll8lli?] and references 
therein. Preheating as a source of gravitational 
waves was first discussed by Khlebnikov and 
Tkachev [20j , and the topic was revived by East- 
her and Lim [S] in 2006, who discussed the scaling 
relationships between the gravitational wave sig- 
nal and the inflationary scale, confirmed in |10| . 
At least three other groups have developed codes 
for computing signal [9I11I12ITB] , which is a non- 
trivial numerical problem, and there is excellent 
agreement between them |13j , although the inter- 
nal workings of the codes differ significantly. Con- 



1 



2 



sider a toy model comprised of a classical inflaton 
(j) coupled to a second field x with Lagrangian 

The equation of motion for x is 



X + ^Hx 



1 



(2) 



Expanding x in terms of its momentum (Fourier) 
modes we obtain 



Xk + 3Hxk +[K+ g^cf' I Xfe = 



(3) 



As inflation ends, cj) oscillates about zero. If we 
ignore the expansion of the universe H = and 
4> oscillates with constant amplitude $, although 
backreaction from x particle creation extracts en- 
ergy from the (j) field, eventually terminating res- 
onance. Changing variables to 



q = -— -, A = — + 2(j, z = mt, 

turns ^ into a Mathieu Equation, 
Xk + iA-2qcos{2z))xk^0. 



(4) 
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where primes denote differentiation with respect 
to z. Every solution to Mathieu's equation has 
two parts, 



Xk oc /(z)e 



(6) 



where f{z) is periodic and /i isthe Mathieu char- 
acteristic exponent. If /i has an imaginary part 
the solution has an exponentially growing mode. 
Figure [T] shows the values of as a function 
of A and q. We see that (|4]) requires A > 2q, so 
we are interested in the parameter values that lie 
to the left of the diagonal line in Figure [1] 

The amplification of selected Xk (and, by back- 
reaction, makes the universe increasingly in- 
homogeneous, yielding an inhomogeneous, time 
dependent energy density - which necessarily 
leads to the emission of gravitational radiation. 
Ultimately, the only way to follow the full evo- 
lution is via numerical simulation, and then ex- 
tracting the contribution to the tensor modes. 




Figure 1. The imaginary part of the Mathieu 
critical exponent is plotted, with darker colors 
corresponding to a larger imaginary component. 
Outside the heavy black lines the exponent is 
real- valued, and the corresponding solutions are 
strictly oscillatory. The diagonal line corresponds 
to A = 2q. 



The stress-energy tensor associated with grav- 
itational radiation is given by |21j 

1 



T, 



(7) 



32ttG 

and is specific to the transverse-traceless part of 
the metric perturbation. The (• ■ •) denotes a spa- 
tial average and since we are solve for the field 
values (and /i^ ) numerically, we simply integrate 
over the full "grid" on which our numerical so- 
lutions are computed. The associated energy en- 
ergy density is the 00 component, 

1 /. .„\ NT^ 1 
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We compute /ly by extracting the source terms 
for the equations of motion obeyed by the from 
the full T^jy obtained from numerical solutions of 
the scalar field dynamics 1 . We express the grav- 
itational wave power in terms of its contribution 
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to the fractional energy density in gravitational 
waves, flgw, per logarithmic interval in wavenum- 
ber, 

dflg-u, _ 1 dp _ nk^ ^^j^ 2 (Q\ 

which we then convert in present-day units via an 
appropriate matching condition [Ij. 

The resulting gravitational wave spectrum is 
sharply peaked, and located at very high frequen- 
cies when compared to both suspected astrophys- 
ical sources, and current and future interferomet- 
ric experiments. The gravitational wave back- 
ground depends on the details of the resonance 
bands in the underlying cosmological model, and 
modes move in and out of resonance bands as the 
universe expands. For the example considered 
here, the lowest lying resonance band is domi- 
nant. Modes whose physical wavenumber k/a{t) 
never enter this region during preheating are not 
significantly amplified, resulting in the cutoff we 
see at high frequencies in Figure [3l 

In this model (and many resonant scenarios), 
the resonant modes are not vastly smaller than 
the post-inflationary Hubble scale. If the uni- 
verse promptly reheats to the GUT scale the co- 
moving post-inflationary Hubble scale is mapped 
to centimeter scales in the present universe, and 
the spectrum peaks near GHz scales. If the in- 
flationary scale is lowered, the signal appears at 
lower frequencies, because the growth of the uni- 
verse between the present day and the preheating 
epoch is necessarily smaller [8] - and if the energy 
scale is very low (closer to the TeV scale than the 
GUT scale), a background of gravitational waves 
generated during preheating would provide an al- 
ternative target for experiments like BBO. 

The amplitude of this signal does not neces- 
sarily change with the inflationary energy scale. 
Only the gradient terms of the scalar field which 
source the gravitational wave background, and 
in many preheating scenarios the gradient energy 
rises until its average value is similar to that of 
the kinetic and potential energies. At fixed den- 
sity this amplitude of the source terms are thus 
maximized. Conversely, the initial power gen- 
erated will decrease with the inflationary scale, 
since the quadruple and higher moments will be 




Figure 2. A cross section of the x fi^ld at three 
times, with an arbitrary normalization. Gravita- 
tional waves are induced by (time varying) spatial 
gradients, which peak in the middle frame. 
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Figure 3. 



The gravitational wave power (at the present epoch) generated with m = 10 
2.5 x 10^. 



and 



smaller. However, the overall expansion of the 
universe between preheating and the present day 
is also reduced in these scenarios. The energy 
density in gravitational waves scales like radia- 
tion, so p)gw ^ {ttpreheat I aoY ~ ^^'^ the Combina- 
tion of these two factors cancels, leading to a max- 
imal amplitude of Vlg^{k) ^ - a number 
also seen in bubble collision calculations [22l23j . 



3. PRIMORDIAL BLACK HOLES 

A second, and more recently explored, mecha- 
nism which may generate gravitational waves in 
the post-inflationary universe is Hawking radia- 
tion from decaying Primordial Black Holes [PBH] 
[2]. This analysis focusses on very small black 
holes, which decay before the onset of nucleosyn- 
thesis, and whose initial abundance is otherwise 
very weakly constrained. It has long been rec- 
ognized that gravitons would be emitted as these 
black holes decay, although the precise rate is sig- 
nificantly modified by grey body corrections [24] . 

Assume a PBH population whose mass is equal 
to the energy contained inside the Hubble volume 
at the instant they collapse. Recalling that = 



8TTp/3Mp, and defining p — E^^^.^, 



M 



BH 



3 M-^ 
327r £;.2 . 



(10) 



which is the mass contained inside a sphere of ra- 
dius \/H. Including grey body corrections Tsi a 
Schwarzschild black hole emits (massless) parti- 
cles with momentum k, reducing its total energy 
as 

(11) 



dE 



2.9 M 



2 

BH 



dtdk 



T 



TT M4 e^/T - 1 

M2 



(12) 

where g is the effective number of (bosonic) de- 
grees of freedom, after integrating over grey body 
corrections [2]. Here k is the physical wavenum- 
ber and the black hole lifetime is 



10, 2407r M| 



BH 
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3 M5 



240 



(13) 



An upper bound on E^^^^ comes from the infla- 
tionary energy scale, which is constrained by the 
non-detection of a primordial gravitational wave 
background in the CMB. Since we are interested 
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Figure 4. ngw{f)h? with (from left to right) 
E,„,, = 10l^ IQi"*, 1013, and lO^^ GeV. In all cases 
13 = 0.001 and g = 1000. 



in black holes which decay prior to nucleosynthe- 
sis with time to spare for thermalization, we need 
T < 100 s. With = 1012 GeV, r « 29/g s, 
and the initial temperature is 18.8 TeV. Standard 
model states alone give g ^ 0(10^). 

The mass- fraction of PBH is denoted fl bh- Ini- 
tially ^BH = /3. We assume /? <C 1 and that the 
remaining matter consists of radiation. Denoting 
the number density of PBH by n{t), the energy 
density psH = n{t)MBH{t) we get the following 
set of equations 



dt 



dt 

a 



= n{t)MBH + 


n{t)MBH , 


(14) 
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along with equation ()11|) from which one obtains 
flgw{k) by an appropriate rescaling. 

We again express our results in terms of the 
present-day spectral energy density of gravita- 
tional radiation [8113] . Figure 2] shows the 
present-day power in gravitational radiation as a 
function of E^^^ . The gravitational wave power is 
substantial, and found at very high frequencies, 



relative both to any preheating signal, and the 
electromagnetic radiation in the CMB. Roughly 
speaking, the temperature of the universe scales 
as l/a{t). A decaying black hole is much hot- 
ter than the surrounding universe, but the emit- 
ted gravitational waves are redshifted by the 
same factor as other radiation. Consequently, 
these gravitational waves will necessarily have 
a higher frequency than the present-day CMB, 
which peaks around 100 GHz. Lowering E^^^ in- 
creases the PBH lifetime, enhancing this discrep- 
ancy and pushing the gravitational wave signal to 
higher (present day) frequencies. The signal also 
depends on g and /3, as explained in [2]. 

These primordial black holes survive for many 
Hubble times, and since their density scales like 
matter, even a small initial population can lead 
to a universe that undergoes a transient period of 
matter domination before the black holes decay, 
which then rethermalizes the universe. The addi- 
tional phenomenology arising from any such mat- 
ter dominated stage has not been fully explored, 
and remains a topic for future enquiry. 

4. DISCUSSION 

I have reviewed two different mechanisms for 
generating gravitational waves at the end of in- 
flation - via parametric resonance or preheating, 
and graviton emission from decaying primordial 
black holes. Physically these mechanisms are very 
different ~ the first is driven by the higher mo- 
ments of the matter distribution and thus op- 
erates in the purely classical limit. Conversely, 
Hawking radiation is intrinsically quantum me- 
chanical, and would not take place if ?i = 0. 
However, both mechanisms depend critically on 
the detailed behavior of the universe makes the 
transition between inflationary and "regular" ex- 
pansion. Parametric resonance and preheating 
depends sensitively on the couplings between the 
inflaton and other particles, whereas primordial 
black holes only form in models with a sharp rise 
in the primordial perturbation spectrum as infla- 
tion comes to an end. 

In most cases, the frequencies of these gravita- 
tional wave backgrounds far exceed those which 
modern experiments such as LIGO or VIRGO, 
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and space-based successors such as LISA, are de- 
signed to detect. Consequently, if meaningful 
constraints are to be placed on the backgrounds, 
new detector technologies are likely to be re- 
quired. The one exception to this statement is 
if preheating occurs after low scale (ie TeV) infla- 
tion, in which case the gravitational wave signal 
could be visible in proposed space-based interfer- 
ometric detectors. This is particularly relevant to 
BBO, which has been designed to detect the pri- 
mordial gravitational wave signal generated dur- 
ing GUT scale inflation. This signal is invisible in 
low scale models, and the preheating signal pro- 
vides an alternative source for a stochastic back- 
ground, widening the set of inflationary models a 
BBO style experiment would constrain. 

Finally, it is clear that these signals are all "fu- 
turistic" in that placing meaningful constraints 
on their amplitude would require second or third 
generation gravitational wave experiments, and 
at the time of writing no direct gravitational wave 
detections have been reported. However, from a 
philosophical standpoint they serve to underline 
the rich phenomenology that can be associated 
with inflationary models, and the importance of 
fully exploring their consequences and possible 
observable fingerprints. 
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